New probe of modified gravity 
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We suggest a new efficient way to constrain a certain class of large scale modifications of gravity. We show 
that the scale-free relation between density and size of Dark Matter halos, predicted within the ACDM model 
with Newtonian gravity, gets modified in a wide class of theories of modified gravity. 



Models with the large scale modification of gravity are ac- 
tively discussed in the recent years in connection with the 
observed accelerated expansions of the Universe and because 
they can be related to the existence of extra dimensions 01 [2]. 
However, the general principles of gauge invariance and uni- 
tarity strongly constrain possible theories of gravity, modify- 
ing the Newton's law at large distances (see analysis of 0]). 
Thus, in addition to its phenomenological applications, this 
problem is related to the fundamental questions of particle 
physics, field theory and gravity. It is therefore important to 
search for large scale modifications of gravity experimentally. 

A possible set of consistent (as a spin-2 field theory) large 
scale modifications of gravity is described by two parame- 
ters - scale r c and a number < a < 1 IJ-Ql]- r c marks 
the distances at which at the linearized level gravitational law 
changes from 1/r 2 to some other power l/r n , and the pa- 
rameter a determines the value of n. Phenomenologically, 
deviations from Newton's law we are looking for may be rep- 
resented in this parameter space. A significant fraction of this 
space is excluded by precision measurements of the Moon or- 
bit @]. Other natural probes of such modifications are cosmo- 
logical observables (see e.g. and refs. therein). 

In this work we identify a new observable sensitive to the 
large scale modifications of gravity. We demonstrate that uni- 
versal properties of individual dark matter halos are also af- 
fected by the modifications of gravity, and provide novel way 
to probe them. Namely, we show that the scale-invariant rela- 
tion between density and size of dark matter halos, predicted 
by Newtonian gravity within the ACDM model lfl3ll and found 
to hold to a good precision in observed dark matter halos lfl4ll . 
may receive non-universal (size-dependent) corrections for a 
wide range of parameters r c and a. 

Formation of structures in the Universe is an interplay 
between gravitational (Jeans) instability and overall Freed- 
man expansion. The gravitational collapse does not start un- 
til the potential energy U of a gravitating dark matter sys- 
tem overpowers the kinetic energy of the Hubble expansion 
JC ~ ^H 2 R 2 . Once the gravitational collapse has began, at 
any moment of time t a dark halo is confined within a sphere 
of zero velocity or a turn-around sphere. As Hubble expan- 
sion rate H(t) decreases with time, the turn-around radius 
i? ta (i) grows. In the Newtonian cosmology with potential 
<t>N{r) = —GM/r the turn-around radius i? ta is 
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(today for masses ~ 1O 12 M the turn-around radius is 



IMpc). Notice that at any moment of time the average density 
within a turn-around radius ((TJ is proportional to the cosmo- 
logical density and is the same for halos of all masses: 
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It was shown in IU3I1 that the property (|2]i leads to a uni- 
versal relation between characteristic scales and densities of 
dark matter halos. This relation holds in wide class of dark- 
matter dominated o bjec t s (f rom dwarf galaxies to galaxy clus- 
ters) lfl4ll (see also lll5l - [l7lo . The relation is in a very good 
agreement with pure dark matter simulations lfl8i fl9fl . sug- 
gesting that baryonic feedback can be neglected in this case. 
Therefore, this relation can serve as a new tool of probing 
properties of dark matter and gravity at large scales. 

The relation (|2]i continues to hold in the Universe where 
gravity is modified by the cosmological constant A. The grav- 
itational energy of a body of mass M at distance r becomes 
U\ = — — ^jr- . Comparing it with the kinetic energy of 
the Hubble flow K, one arrives once again to the relation (|2]i 
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(c.f. IU3I1 ). The relation (Q~|i still holds and is again independent 
on the mass of the halo. 

What is the most general form of gravitational potential, for 
which the property (0 remains true? Clearly, it will hold for 
all the gravitational potentials of the form 
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where is some constant with dimension of density. In par- 
ticular, A-term obeys this property (with p± oc A/G). All the 
theories of the form (O obey the property that relative cor- 
rection to the Newtonian potential 0jy depends only on the 
density p(r) within a radius r (and not on the mass or the size 
of objects). 

Next, we consider the modifications of gravity 00]. The 
gravitational potential of a spherically symmetric system of 
mass M there has the form 
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where < a < 1. Here the characteristic (Vainstein) radius 
rv is defined as 0,11 |U 
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where f3 = 4(1 — a) 
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If the scale r c is of the order of ~ Hq , such modifications 
of gravity can provide an explanation for the late-time cosmo- 
logical expansion of the Universe lHIH]- The corrections to 
Newton's law become negligible as r — > (7r(0) = 1) and the 
radius (O characterizes the scale where the deviations from 
Newton's potential become of order unity. Using the relation 
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(where p = Af/r 3 ) we find that among the theories of mod- 
ified gravity © only p = 2 (a = \, the DGP model £0]) 
possess the property (|4]i and consequently (0. 

The property ^ can be probed experimentally. Exten- 
sive catalog of DM-dominated objects of all scales, collected 
in (3, exhibits a simple scaling relation of the properties of 
the DM halos. Dark matter distribution in the majority of ob- 
served objects can be described by one of the universal DM 
profiles (such as e.g. NFW JH or Burkert JH). Such pro- 
files may be parametrized by two numbers, directly related 
to observations - a characteristic radius rc (off-center dis- 
tance where the rotation curve becomes approximately flat, 
equal e.g. to r s for NFW) and a DM central mass density 
pc, averaged inside a ball with the size rg. It was shown 
that DM column density S oc pc r c, (see GHHH for a de- 
tailed definition) changes with the mass as S oc M K , where 
K ~ 0.22 — 0.33. It was demonstrated in |13] that in the 
simplest self-similar model (i.e. assuming that rc/Rt& is the 
same for the DM halos of all masses) property (QJ[2]i implies 
a scaling S oc M 1 / 3 , compatible with observations. Observa- 
tions (see e.g. j23)-^3]) demonstrate that the ratio of rc to the 
virial radius depends weakly (as M"" 0,1 ) on the mass of DM 
halos. The Fig. Q] shows ratio of the virial radius of a halo to 
its rc for DM density profiles from the catalog of 11411 . These 
results are in perfect agreement with the ACDM numerical 
simulations (see e.g. Owl l28tn . Due to this slight deviation 
from the self-similarity, the best fit value of the scaling pa- 
rameter k = 0.23 (see IU3I1 for discussion). For the qualitative 
discussion of this work, it is important that S(M) is a fea- 
tureless power-law dependence, whose slope does not depend 
on mass and that the deviation from the slope | is small (as 
follows from observations). 

We can conclude that in the theories, satisfying the con- 
dition (|5]l (e.g. in DGP model) the properties of DM halos 
theory will follow the same scaling relation S oc A/ 1 / 3 and 
the difference with the ACDM case will only be in a different 
normalization of this scaling relation. For other theories, de- 
scribed in |3|,|4j], with a ^ i, we can see from the Eq. (0) that 
the potential <j> a {r) is not of the form and we can expect 
deviation from the universal scaling law. 



I. S — M RELATION FOR GENERAL a 

Let us work out the S — M for a general a in details. A 
general expression, relating the turn-around time, turn-around 
radius and mass within this radius follows from energy con- 
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FIG. 1: Comparison of c v i r = Rvii/rc as a function of DM halo 
mass for observed DM density profiles from the catalog of Jl4ll . 
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Using the general form (0 we can rewrite the expression ([8]) 
in the following form: 
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where dimensionless ratio x ta = Rt&/rv and the function 
I(x ta ) is given by 
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The solution of this equation gives us the "density" p ta = 
A//i? 3 a as a function of M, When r c — s- oo the turn-around 
density p ta becomes 
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Here the constant po is a function of lifetime of the Universe 
only and does not depend on parameters of a dark matter halo. 
This gives a desired relation between a turn-around density 
and the life-time of the Universe in the pure Newtonian cos- 
mology (without cosmological constant). The function I(x) 
is defined in such a way that in the Newtonian limit ir(x) = 1 
one gets I(x) = 1. 

The derivation of Eq. (0 demonstrates that for all theories 
of gravity of the form (0]i (including A-term and the DGP 
model) p ta is a function of po only and does not depend on 
the mass/size of a particular halo. As a result S oc M 1 / 3 
(see HI] for details). 

Further analysis of Eq. © depends on the form of the func- 
tion tt(x). Its exact form is not known (apart from the DGP 
case). Let us start with analyzing several limiting cases. 
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FIG. 2: Function Ii(a) defined in Eq. J 1 3 b 
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If the Vainstein radius ry <C R ta for halos of all masses 
that are experimentally observed (roughly from ~ 10 M 
to ~ 10 16 Mq), then for distances r ^> ry the corrections 
to the Newtonian potential reduce either to the order one 
renormalization of the gravitational constant (on the "normal 
branch") or become indistinguishable from the A-term ("self- 
accelerated branch"). In both cases S(M) oc M 1/3 with the 
normalization, different from the pure Newtonian case. 

In the opposite case ry 3> Rta, one can utilize the per- 
turbative expansion of the function tt(x). The gravitational 
potential well inside the Vainstein radius is given by Jl,01 
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FIG. 3: Examples of <S — M relation in theories with ot < 1/2 
together with the data from fl4ll . 



Under the assumption of exact self-similarity, discussed 
above (i.e. rc/Rta = const) one arrives to the following 
expression for S (recall that j3 = 4(1 — a)): 
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where 



where ci ~ 0(1) and is positive for the "normal branch" and 
negative for the "self-accelerating branch" Jit] in full analogy 
with the DGP model. Notice that a > 1 for a < § . Using 
expansion (fT2l one arrives to 
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where the function I±(a) is shown on the Fig. [2] 

Substituting the expression ( fT3l ) back into equation ^ and 
using (O, we obtain 
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As Xta "C 1 and a > 1, one finds that 



Pta ~ p (l + ci7i(a).< a (p )) 
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where to compute x ta we use Eq. (O with p Q instead of p ta . 
From Eqs. (Q and (TT~5T > we see once again that for all a ^ h 
(i.e. (3 7^ 2) the turn-around density p ta loses its univer- 
sality and becomes the function of the halo mass M. The 
turn-around radius R ta (M) is related to p ta via i? ta (M) = 
(M/ Pta ) 1 /3. 
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An example of the relation ( fTTT i for several a's and r c is shown 
in Fig. |3] 

Clearly, the most interesting case is when ry ss R ta for 
some range of observed halo masses. In this regime the de- 
viations from Newtonian gravity become the strongest. The 
range of values r c for which this happens is shown in Fig. H] 
(the value of to is chosen to be the lifetime of the Universe 
to ~ 1.3 x 10 10 years). We expect that for r c in the region 
Fig|4]the slope of the relation S oc M K will change. Anal- 
ysis of this case requires however an exact solution of the 
non-linear analog of the Poisson equation in theories with a 
(see e.g. 0]), i.e. the knowledge of properties of the func- 
tion it (r fry) in the range of radii, where the perturbative 
expansion (fT2l breaks down. Notice, that the region where 
ry w i? ta shrinks toward the value r c = -to as a |. For 
this value of r c the Vainstein radius in the DGP model is equal 
to the turn-around radius for all halo masses. 1 Knowledge of 
the potential tt(x) would also allow to probe the modifications 
of gravity directly in the Local Group and other nearby galax- 
ies, by studyin g th e infall trajectory around the turn-around 
radius (see e.g. 12911 '). 



In general for r c ~ H 1 the turn-around radius is smaller than ry for 
a < | and bigger than ry for a > | . 
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FIG. 4: Range of r c for which vy ~ i?ta for some halo masses in 
the range 10 10 - 1O 15 M . 

Another way to probe the S — M relation for general r c and 
a is to do numerical simulations in the theories of modified 
gravity (see examples in lflll[l2lo and compare directly with 
observations both the scaling of the central values of S and M 
and the scatter around it. 

Conclusion. The main purpose of this work was to iden- 
tify a new observable that can be used to constrain the large 
scale modifications of gravity. We see that the scaling prop- 
erties of dark matter halos are sensitive to such modifications. 



We demonstrated that the models with predict the de- 

viation from a simple power-law scaling in the S(M) rela- 
tion. Comparison of predictions of such models with the data, 
collected in jTill potentially allows to restrict the values of 
r c from below for a given a. The improved data-processing 
and new observational data on DM distributions will allow to 
strengthen these bounds and make them quantitative. 

In this work we have analyzed only the case when the turn- 
around sphere is well inside the Vainstein radius, ry ^> i? ta . 
To analyze a general case, a better theoretical understanding 
of the function ir(r) (defined via Eq. ©) is needed. Together 
with better quality of data this will allow to extend our analy- 
sis to a wider range of parameters. 

In the case a = \ (the DGP model) the S(M) dependence 
remains featureless. In this case one has ry sa R ta for all 
masses (for r c ~ H^ 1 ) and the deviations from the Newto- 
nian gravity at turn-around radius will be strong. Therefore 
this model (in general, all models that have ry ~ R ta for ha- 
los of M ~ 10 12 M Q ) can be probed by studying the infall 
trajectories around the turn-around radius in the Local Group 
and nearby galaxies l29ll using the available data and the data 
from forthcoming surveys of the Milky Way as well as GAIA 
mission. 
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